Electronic correlations and crystal structure distortions in BaBiC>3 
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BaBiC>3 is a material where formally Bi 4+ ions with the half-filled 6s-states form the alternating set 
of Bi 3+ and Bi 5+ ions resulting in a charge ordered insulator. The charge ordering is accompanied 
by the breathing distortion of the Bi06 octahedra (extension and contraction of the Bi-O bond 
lengths). Standard Density Functional Theory (DFT) calculations fail to obtain the crystal structure 
instability caused by the pure breathing distortions. Combining effects of the breathing distortions 
and tilting of the BiC>6 octahedra allows DFT to reproduce qualitatively experimentally observed 
insulator with monoclinic crystal structure but gives strongly underestimate breathing distortion 
parameter and energy gap values. In the present work we reexamine the BaBiC>3 problem within 
the GGA+U method using a Wannier functions basis set for the Bi 6s-band. Due to high oxidation 
state of bismuth in this material the Bi 6s-symmetry Wannier function is predominantly extended 
spatially on surrounding oxygen ions and hence differs strongly from a pure atomic 6s-orbital. That 
is in sharp contrast to transition metal oxides (with exclusion of high oxidation state compounds) 
where the major part a of d-band Wannier function is concentrated on metal ion and a pure atomic 
rf-orbital can serve as a good approximation. The GGA+U calculation results agree well with 
experimental data, in particular with experimental crystal structure parameters and energy gap 
values. Moreover, the GGA+U method allows one to reproduce the crystal structure instability due 
to the pure breathing distortions without octahedra tilting. 

PACS numbers: 



I. INTRODUCTION 



BaBiOa is a parent compound for high-T c supercon- 
ductors Bai-ajK^BiOa and Bai-^Pb^BiOs. Properties 
of these compounds are the most intriguing when chem- 
ical composition is near transition to SC state. It would 
be promising to describe softening of phonon modes 
with doping as a result of ab initio calculation in order 
to support phonon mechanism for superconductivity in 
these materials. The proper description of electronic and 
crystal structure of pure BaBiC>3 has serious difficulties 
within Density Functional Theory (DFT). A possible so- 
lution to this problem is to take into account strong elec- 
tronic correlations that are important in high-T c cuprates 
and oxipnictides. However, the most straightforward way 
to include correlation effects: the LDA+U method fails 
in attempts to improve results of the DFT calculation 
for BaBiC>3. In the present work it is shown that the 
reason for this failure is a high oxidation state of the 
compound. A conventional basis set used in the LDA+U 
method is pure atomic orbitals, that is usually a good 
approximation for d-bands Wannier functions in transi- 
tion metal compounds (with important exception of high 
oxidation state materials). For BaBiC>3 one should use 
genuine Wannier functions as a basis set for the LDA+U 
correction potential, because they differ a lot from pure 
atomic Bi 6s-orbitals. Below it is demonstrated that the 
LDA+U correction in Wannier function basis results in 
a strong improvement of the calculated parameters for 
crystal structure and spectral properties for BaBi03. 

The LDA+Ul and LDA+DMFT- methods for the 



electronic structure calculations of strongly correlated 
systems are designed as a combination of DFT and model 
Hamiltonian approach (such as Hubbard or Anderson 
models). A DFT potential is defined by electronic den- 
sity only and hence is invariant under unitary transforma- 
tion of auxiliary Kohn-Sham orbitals (DFT knows noth- 
ing about basis functions of a particular method; DFT 
knows only Kohn-Sham orbitals and total charge den- 
sity). However, in model Hamiltonians Coulomb interac- 
tion between electrons is explicitly defined for localized 
states having symmetry of atomic orbitals. A specific 
form of these orbitals remains to be determined in a cal- 
culation scheme developed for LDA+U method realiza- 
tion. For the LMTO^ method, the basis set (Muffin-tin 
orbitals) is explicitly constructed in the form of atomic- 
like orbitals and so in this case the choice of orbitals for 
Coulomb interaction term is obvious. For the LAPW~ 
method, there are atomic symmetry wave functions in- 
side muffin-tin spheres and there are pseudoatomic or- 
bitals even for pseudopotential methods. These orbitals 
are usually used to define the Coulomb interaction term 
in a LDA+U calculation scheme based on corresponding 
methods. 

The most general and "rigorous" way to define atomic- 
like orbitals in solid is to build Wannier functions. That 
is done via Fourier transformation from dclocalizcd Bloch 
functions V'k in reciprocal space of wave vectors k to func- 
tions W T in direct space lattice sites vectors T. Both 
sets, the Bloch functions -0k and Wannier functions W T , 
generate the same wave functions Hilbert space and can 
be considered as different basis sets for this space. Cor- 
relation effects are usually strong enough to be taken 
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into account for transition metal compounds, where d- or 
/-states produce narrow energy bands. Fourier transfor- 
mation for corresponding Bloch functions results in Wan- 
nier functions, which are well localized inside an atomic 
sphere and hence using atomic d- or /-orbitals instead of 
them is a good approximation. 

Correlated states differ substantially from pure atomic 
d-orbitals in compounds with elements in high oxidation 
state, for example, CaFe03, AuTe2, and BaBi03. Nar- 
row energy bands near the Fermi level have a serious 
admixture of p-states in such compounds, which could 
be even larger in value than d-states contribution. More- 
over, an energy gap could be formed predominantly by 
p-states with a slight admixture of d-states (as it happens 
in AuTe2^. The mixed nature of correlated electronic 
states does not allow to use traditional implementations 
of the LDA+U method with pure atomic orbitals as a ba- 
sis set for Coulomb interaction term. Explicit Wannicr 
functions basis set for correlated states should be cho- 
sen to perform electronic structure calculations for these 
materials. 

Many researchers used Wannier function basis for 
strongly correlated materials calculations in recent years, 
for example, see Refs^— . In the present work we inves- 
tigate BaBiOa where Wannier functions basis in contrast 
to atomic orbitals is crucial for proper crystal structure 
description. In the case of BaBiC>3 partially filled band, 
which is a source of correlation effects, is formed by Bi 
6s-states and oxygen 2p-states. The energy of original 
Bi 6s-level is below oxygen band but strong hybridiza- 
tion between them results in the splitting of the partially 
filled band upward from the 2p-states manifold. The con- 
tribution of oxygen states to this band is stronger then 
contents of Bi 6s-states but corresponding wave function 
still has s-symmetry in respect to Bi site. Wannier func- 
tion constructed from this band has the same symmetry 
but it is spatially distributed more on oxygen sites than 
on Bi ion. Accounting for the Coulomb interaction on 
such a Wannier functions basis leads to the results very 
different from the LDA+U calculations using pure Bi 6s- 
orbitals. 



II. CALCULATION METHOD 

A. Wannier functions 

Wannicr functions (WFs) |Wjf ) arc defined as Fourier 
transforms of the Bloch functions | ^nk) : 

i^ T > = ^E e " 4kT i^>, a) 

where T is the lattice translation vector, n is the band 
number, and k is the wave vector. WFs seem to be 
an optimal choice to describe correlated states because 
these functions are localized, site-centered and repre- 
sent a complete basis set for the Bloch functions Hilbcrt 



space. In the present paper WFs were generated as pro- 
jections of Bloch sums of the atomic orbitals \4> n k.) = 
^ T e lkT |^) onto a subspacc of the Bloch functions 
(the detailed description of WFs construction procedure 
within pseudopotential method is provided in RefJ^): 




N 2 

\W n ±) = 2 |*<k>(*<k|^k>. (3) 

f=JVi 

The resulting WFs have symmetry of the atomic or- 
bitals 4> n and describe electronic states that form energy 
bands numbered from N\ to N-2- To describe Coulomb 
repulsion between occupied and empty states one needs 
to compute occupation numbers of WFs. The Wannier 
functions occupancy matrix Q^m ^ s given by 

Qui = M [Y, E - E i)(**\ J K) 

V k i=Ni / 

(4) 

where 9 is the step function, £j(k) is the one-electron 
energy for the state i, and Ef is the Fermi energy. 

B. LDA+U method in Wannier functions basis 

One of the most general and accurate approaches to de- 
scribe the electronic structure of strongly correlated sys- 
tems is the dynamical mean-field theory (DMFT)^i~— . 
The method combining ab initio DFT approach and 
model DMFT calculations - LDA+DMFT was applied 
recently to the quantitative description of both electronic 
and crystal structure of several compounds with strong 
electronic correlations^^. The LDA+DMFT method 
is successful and promising but it demands a lot of com- 
putational resources. In the present paper the static limit 
of Dynamical Mean-Field Theory - the LDA+U method 
is used. There is an additional term in Hamiltonian op- 
erator to take into account Coulomb correlations. If one 
considers on-site interaction only and takes into account 
the diagonal form of the occupation matrix Qj^f , the 
additional term is written in WFs basis as 

ff'(k) = \W mk ) (u*(\- fj) (Wmk\, (5) 

where n m = Q^m ~ the occupation number for m-th WF. 
In the case of BaBi03 the non-diagonal terms of QjJ^f are 
equal to zero by construction (we have strictly orthogonal 
WFs of Bi-s-orbital symmetry) . For more complex basis 
sets it could be necessary to use rotationally invariant 
potential correction and take into account non-diagonal 
terms. 

The full Hamiltonian operator is written as 

H LDA + U = H LDA + H', (6) 



where H LDA - the Hamiltonian operator within LDA (or 
GGA) approximation. The form of the correction poten- 
tial ((5]) results in a negative addition to potential —U/2 
for occupied states n m = 1 and in a positive value +U/2 
for empty states n m = 0. Hence the LDA+U correction 
to DFT will increase energy separation between bands 
below and above the Fermi level increasing energy gap 
values and enhancing total energy gain of the distortion 
that causes the gap appearance. 

The value of Coulomb interaction parameter U is com- 
puted via constrained LDA calculation as described mi£. 

The total energy is computed as 

E tot = E LDA + E u - E DC , (7) 

where E LDA is the total energy from a standard DFT 
calculation (LDA or GGA could be used, we preferred 
GGA), E u = lE m ^ m/ Un m n m /, E DC = \Un(n - 1) is 
the double counting correction, and n = Tr(Q^ AV J^) is the 
total occupancy of WFs. 



III. RESULTS AND DISCUSSION 

Undopcd BaBiOa has a monoclinic crystal structure 
(symmetry group is C2/m). The structure could be ob- 
tained from an ideal cubic perovskite structure simul- 
taneous breathing distortion of the BiOg octahedra and 
tilting the octahedra around the [110] axis. Primitive cell 
contains two formula units. Chemical formula could be 
written^ as Ba^ + Bi 3+ Bi 5+ 0^". Two different lengths of 
Bi-0 bonds correspond to two different valences of Bi ion: 
the short to Bi 5+ -0 and the long to Bi 3+ -0. The neigh- 
boring Bi-type ions for Bi 3+ are Bi 5+ and vice versa. The 
alternation of Bi 3+ and Bi 5+ ions forms charge density 
wave. Experimental valued for the breathing distortion 
is 6=0.085 A and for the tilting is t=10.3°. Spectroscopy 
measurements show that BaBi03 is an insulator with the 
energy gap value «0.5 eV— . 

Liechtenstein et al. performed calculations of elec- 
tronic and crystal structure of BaBi03 with full-potential 
LMTO method^ . The calculated values of the Bi0 6 
octahedra distortion (6=0.055 A, t=8.5°) underestimate 
the experimental data as well as the energy gap value 
«0.02 eV. Also authors were unable to obtain cubic crys- 
tal structure instability against pure breathing octahedra 
distortion without tilting. 

Franchini et al. have modeled the crystal structure dis- 
tortions more successfull y 21 ' 22 . The authors used Heyd- 
Scuseria-Ernzerhof (HSE) hybrid functional that includes 
25% of exact Hartree-Fock exchange in addition to DFT, 
so it partially takes into account correlation effects. The 
hybrid functional allowed to reproduce energy gap value 
with slight overestimation (0.65 eV) and get octahedra 
distortion values close to experimental ones (6=0.09 A, 
i=11.9°). 

Thonhauscr and Rabe— have reported the cubic crys- 
tal structure instability in respect to breathing distor- 
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FIG. 1: (color online) Squared moduli of the Wannier function 
centered on Bi 3+ ion (a big red sphere in the center). The blue 
spheres are oxygen ions. 



tion of the BiOg octahedra using obsolete Local Den- 
sity Approximation ( LDA) while more rigorous Gener- 
alized Gradient Approximation (GGA) calculations do 
not show any minimum at all (see Fig. [3]). The total 
energy lowering value for distorted structure in LDA cal- 
culations— was 3.35 meV per formula unit corresponding 
to crystal structure transformation temperature as low as 
30 K, that clearly contradicts to experimental fact that 
distorted cubic crystal structure of BaBiOa is stable at 
room temperature and higher temperatures^. 

In the present work to describe the correlated states 
with WFs, a self-consistent calculation within the gen- 
eralized gradient approximation (GGA) was performed 
as the first step. The pseudopotential plane-waves 
method implemented in the Quantum-ESPRESSO pack- 
age^, was used. Vandcrbilt ultrasoft pseudopotentials^ 
were taken from QE pseudopotentials library. A kinetic- 
energy cutoff for the plane-wave expansion of the elec- 
tronic states was set to 45 Ry. Integrations in reciprocal 
space were performed using (10,10,10) Monkhorst-PackSi 
k-point grid in the full Brillouin zone. 

The GGA calculation results in semiconducting solu- 
tion with a very small indirect energy gap «0.05 eV. Two 
energy bands in the interval [-1.5; 2.5] eV, see FigJH are 
generally formed by the O-p states, but they have symme- 
try of the Bi-6s states. The Bi-5d and Bi-6p partial den- 
sities of states give vanishing contribution to the states 
near the Fermi level. 

WFs were generated by projection of two Bi-s atomic 
orbitals onto subspace defined by two energy bands near 
the Fermi level. Squared moduli of WF centered on Bi 3+ 
ion is shown in Fig. [TJ WF has symmetry of s-orbital. 
However, it has significant contribution from p-states of 
neighboring oxygen ions. The second WF (centered on 
Bi 5 + ion) will differ only by the ratio between O-p and 
Bi-s states. 

The contribution of oxygen states to partially filled 
bands could be accounted for by including in the Hamil- 
tonian Hubbard U correction to the oxygen p-states in 
addition to Bi-s orbitals. However in that case one needs 



4 



to take into account also off-diagonal intcr-atomic s — p 
terms in occupation matrix. The WF approach proposed 
in the present work is simpler. 

To perform calculation within LDA+U method in WFs 
basis, the U value for BaBi03 was computed via the con- 
strained LDA method. The Hubbard U value could be 
calculated as the second derivative of total energy over 
occupation of correlated stated: 



U = 



J DFT 



d 2 n c 



(8) 



In DFT one-electron energy of a state is a derivative of 
total energy in respect to the state occupancy. Then U 
value can be calculated as: 



U 



de c 



(9) 



where e corr is one particle energy of correlated state, and 
n C orr — its occupation number. The partial derivative 
is taken numerically as described below. Constrained 
potential that affects only states considered as correlated 
ones is defined as a projection operator: 



(10) 



where 6V n is a small number, for example, ±0.1 eV, 
positive for WF centered on one Bi-site and negative for 
WF centered on another Bi-site. The constrained po- 
tential is included into the Hamiltonian operator dur- 
ing self-consistency cycle, that allows to take into ac- 
count screening of charge redistribution. As a result, the 
WFs occupation numbers and energies change from un- 
constrained results. The derivative © is calculated as a 
ratio of these variations. We have checked that the cal- 
culated derivative value of U does not depend on 5V n , by 
using its various values equal to (0.1, 0.2, 0.3, 0.4 eV ). 

The resulting U value for BaBiC>3 equals to 0.7 eV. The 
value of Coulomb parameter strongly depends on spacial 
form of basis function^. Therefore the rather small value 
0.7 eV obtained in our calculations is reasonable for WFs 
shown on Fig. [T] that is rather extended in space. The 
atomic Bi-s orbital is much more localized in space than 
WFs hence the U value for atomic orbitals was found to 
be significantly larger U= 2.5 eV. 

The U value for WFs basis is twice smaller than the 
width of the partially filled energy bands that indicates 
relative weakness of Coulomb interaction in BaBi03. 
But, as it is shown below, this Coulomb interaction plays 
the important role in the formation of the electronic 
states and crystal structure of the compound. 

A simulation of the oxygen octahedra tilting and 
breathing distortion was performed in a monoclinic 
(I2/m) supercell with 2 formula units. Sixty primitive 
cells of Ba2 + Bi 3+ Bi 5+ Og~ were generated correspond- 
ing to different combinations of the breathing and tilting 
distortions. For every structure a self-consistent calcula- 
tion within CCA was done. Then WFs were constructed 
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FIG. 2: Band structure of BaBiOa obtained for various crys- 
tal structures within the GGA and GGA+U approaches, 
(a) Ideal cubic crystal structure, GGA. (b) Distorted crys- 
tal structure with the breathing-only distortions, 6=0.075 A, 
GGA. (c) Distorted crystal structure with the breathing-only 
distortions, 6=0.075 A, GGA+U in WFs basis, (d) Mono- 
clinic crystal structure (6=0.075 A, t=12°), GGA. (e) Mono- 
clinic crystal structure (6=0.075 A, £=12°), GGA+U in WFs 
basis. Zero energy corresponds to the Fermi level. 



and self-consistent calculations within LDA+U (actually 
GGA+U) were performed. Coulomb interaction correc- 
tion was applied only to the states described with the 
two WFs as defined in Eq. (|SJ). 

In Fig. [2] the formation of the energy gap due to the 
cubic structure distortion is shown. Band structure of 
ideal cubic BaBiOa, see Fig. [5] (a), is definitely metallic. 
Two energy bands in the [-1.5; 2.5] eV energy interval 
cross the Fermi level. These bands are double degenerate 
along the A— >-M direction. The cubic crystal structure is 
unstable against distortion that lowers the symmetry and 
opens an energy gap in band structure. 

In BaBiOs the breathing distortion of the oxygen oc- 
tahedra plays the crucial role. In Figs.LH(b) and (c) band 
structure is shown for the crystal with the breathing-only 
distortions of the oxygen octahedra. The breathing value 
is 6=0.075 A. The degeneracy breaking along the A— s-M 
direction is clearly seen in the GGA calculation result. 
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Breathing distortion (A) 

FIG. 3: (color online) Calculated within LDA (magenta 
dashed), GGA (black solid), atomic orbitals GGA+U (green 
dashed), and GGA+U in WFs basis (red solid) total energy 
of BaBi03 as a function of pure breathing for tilting = 0°. 
The inset shows results of GGA+U in WFs basis calculations 
close to zero breathing distortion. 



Two separate bands appear in the [-1.5; 2.5] eV interval 
but the Fermi level crosses both of them and the energy 
gap docs not open. Accordingly there is no total energy 
minimum due to the breathing distortion in the GGA 
calculation. That is clearly illustrated by the black curve 
in Fig. [3J The dependence of total energy of BaBi03 
cell on the breathing distortion calculated within GGA 
has a minimum for the distortion parameter value 6=0 A 
(stable cubic structure). Hubbard interaction correction 
(GGA+U) in the basis of WFs leads to appearance of 
the band gap equal to 0.3 cV, Fig. [2] (c). Also cubic 
crystal structure of BaBi03 is unstable against the pure 
breathing distortion in the GGA+U calculation. The to- 
tal energy of the cell has a minimum for 6=0.075 A, see 
Fig. [3] The energy lowering is equal to 326 meV per 
formula unit. 

The total energy dependence calculated with the 
GGA+U method in atomic orbitals basis is also shown 
in Fig. [3] with the dashed green curve. The minimum of 
the total energy corresponds to the cubic crystal struc- 
ture. One can see that atomic orbitals basis potential cor- 
rection calculation docs not result in an improvement of 
crystal structure description. The energy dependence ob- 
tained with the LDA calculation is shown on Fig. [3] with 
the dashed magenta curve. It has very shallow minimum 
for non zero breathing distortion in agreement with^, 
but there is no any minimum in the curve obtained in 
GGA that is supposed to be more accurate approxima- 
tion than LDA. 

In the monoclinic cell of BaBiO,3 with simultaneous 
breathing and tilting of BiOg octahedra the energy gap 




Breathing (A) 



FIG. 4: (color online) Total energy of BaBiOg as a function 
of breathing and tilting distortions calculated within GGA 
and GGA+U. Zero energy level corresponds to cubic struc- 
ture (6=0, t=0) within GGA. Experimental minimu m 18 ' 24 is 
located at (6=0.085 A, t=10.3°) and is marked by the cross; 
The GGA energy minimum value is at (6=0.04 A, t=12°) and 
the GGA+U value is at (6=0.075 A, t=12°) 



(equal to 0.05 eV) exists in the result of the GGA calcu- 
lation, see Fig.[2](d). Despite the fact that the gap value 
is very small the cubic structure is unstable against com- 
bined distortion in GGA. For the monoclinic structure 
the total energy minimum corresponds to the distortion 
parameters values (6=0.04 A, t=12°). The values under- 
estimate experimental data in agreement with the previ- 
ous calculations^. 

The GGA+U correction for monoclinic crystal cell in- 
creases the energy gap value to 0.55 eV, see Fig.[2](e), in 
a good agreement with the experimental value of 0.5 eV. 

The splitting of the half-filled double degenerate Bi- 
s band and the opening of the energy gap results in a 
filled band corresponding to the 6s-state of one Bi-sitc 
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(Bi 3+ ) and an empty band formed by the states of an- 
other Bi-site (Bi 5+ ). A charge difference between these 
two Bi-sites (equal to 2e) is expected for fully ionic pic- 
ture. However, the occupancy values difference for two 
atomic Bi-6s orbitals in the LDA+U calculation is equal 
to 0.43e only in agreement withal. The small value of 
the charge difference can be understood from an analysis 
of WFs contents. The Bi-6s symmetry Wannier function 
for Bi 3+ site has only 23% of Bi-s atomic orbital and 77% 
of O-p orbitals; on Bi 5+ site WF composition is 20% of 
Bi-s and 80% of O-p. A change in occupation number of 
WF by ±1 results in a change of s atomic orbital occu- 
pancy by +0.23 for Bi 3+ and -0.20 for Bi 5+ , that gives 
0.43 in total. 

In Fig. H the total energy of BaBi0 3 in GGA and 
GGA+U as a function of breathing and tilting distortions 
is shown. Zero energy value corresponds to the GGA to- 
tal energy of ideal cubic crystal structure. The surface 
was obtained by extrapolation from 60 points as the 4 th 
order polynomial. The experimental values of the distor- 
tions are (6=0.085 A. £=10.3°) and is marked by the cross 
in Fig. 2J The GGA calculations gave total energy mini- 
mum for distortion value parameters (6=0.04 A, t= 12°) 
strongly underestimating the experimental data. The 
total energy minimum within the GGA+U calculation 
corresponds to distortion parameter values (6=0.075 A, 
t=12°) that are very close to the experimental ones. 

Figure [4] clearly illustrates the effectiveness of the cal- 
culation scheme proposed in the present paper. The us- 
age of the WFs basis for the partially filled states in the 
GGA+U method allows one to reproduce not only the 



energy gap value, but also noticeably improve the de- 
scription of the crystal structure distortions in BaBiOs. 
The most important result is the instability of cubic crys- 
tal against the pure breathing distortion obtained in the 
GGA+U calculations. 



IV. CONCLUSION 

The scheme for studying of electronic correlations in 
high oxidation stated compounds based on the LDA+U 
method in Wannier functions basis was proposed and 
applied to the case of BaBi03. The energy gap value 
and crystal structure parameters of monoclinic BaBiOs 
were successfully described in agreement with experimen- 
tal data. Instability of cubic crystal structure in respect 
to pure the breathing distortion was obtained for the first 
time. 
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